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Abstract 

A study for checking validity of the auxiliary field method (AFM) is made in quantum 
mechanical four-fermi models which act as a prototype of models for chiral symmetry 
breaking in Quantum Electrodynamics. It has been shown that AFM, defined by an in- 
^ i sertion of Gaussian identity to path integral formulas and by the loop expansion, becomes 

more accurate when taking higher order terms into account under the bosonic model with 
a quartic coupling in 0- and 1-dimensions as well as the model with a four-fermi interac- 
tion in 0-dimension. The case is also confirmed in terms of two models with the four-fermi 
interaction among N species in 1-dimension (the quantum mechanical four-fermi models): 
higher order corrections lead us toward the exact energy of the ground state. It is found 
that the second model belongs to a WKB-exact class that has no higher order corrections 
other than the lowest correction. Discussions are also made for unreliability on the con- 
tinuous time representation of path integration and for a new model of QED as a suitable 
probe for chiral symmetry breaking. 
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1 Introduction 



Chiral condensation, characterized by that the expectation value of fermi field operators, 
(tpip), becomes nonzero, is a nonperturbative phenomenon: in relativistic field theory a 
model with four-fermi interactions proposed by Nambu-Jona-Lasinio .1. has been served 
as a simple example exhibiting chiral symmetry breaking. The Lagrangian reads 

2 

c = -trM + 9 - {(v^) 2 - (?75^) 2 } • (i) 

Originally the mean field had been adopted to deduce the gap equation 1 , but later a 
more powerful prescription has been dispensed by Gross-Neveu then Kugo and Kikkawa 
[2]: the method is evolved by inserting an identity in terms of the Gaussian integration 
with respect to a fictitious field called an auxiliary field or the Hubbard-Stratonovich field 
[3] into a path integral expression of Eq. lfl|). (As for the path integral expression see jl].) 
The prescription combined with the loop expansion [5] is simple and transparent and is 
called the auxiliary field method (AFM) jSJ. (Any approximation schemes cannot spread 
out widely, unless those are simple and transparent. The Feynman graph technique is a 
typical example.) If the number of fermion N is large enough quantum corrections of the 
auxiliary field is negligible. The gap equation obtained under this assumption exhibits 
chiral symmetry breaking. However if N is small, quantum corrections, that is, higher 
loop contributions may change the phase structure 0. 

A question would, then, be raised that how effective is AFM in quantum mechanical 
systems? Accordingly we have investigated the role of auxiliary fields in the bosonic case 
with a quartic interaction in 0- and 1-dimension [Sj and in 0-dimensional four-fermi case 
[2]. These models can be solved exactly or numerically so that we can check how accurate 
is the result of AFM. We find that in the bosonic case AFM does work excellently when we 
take higher loops into account. (Of course, the loop expansion is the asymptotic expansion 
so that we should stop considering higher loops somewhere.) 

In order to illustrate AFM, an outline of the 0-dimensional fermionic model is given 
the target quantity is the "partition function" Z (although in dimension), 



Z = I d N £d N C exp 



where 

N 



(2) 



i=l 



and the coupling constant A 2 is supposed real. We have introduced 2A r -Grassmann vari- 
ables and the notation is followed from the textbook of ref. |1U| . Introduce an auxiliary 
field, y, to kill the (£* -£) 2 term in Eq. (j2j), which can be realized by inserting the identity 



■ :°° dy 
1=1 ^ 6XP 



— oo 



(4) 
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to give 



-oo V 27T 

oo 



y~_ 

' 2 




dy / A 
= \uj H — =y I exp 

-oo v2^ V V^V 



,21 



dyexp[-iV/(y)] 



(5) 



with 



dl 



J(y) = ^--lm> + Ay) 



(6) 



Find a solution y of — = , called the classical solution or the saddle point, then expand 
dy 

I(y) around yo 

r(2) A3) 

m = io + - vof + -^-(y - yo) 3 + ■■■ ; 4 n) = i (n) (yo) ; (7) 



and finally make a change of variable such that (y — yo) 



N 



to obtain 



Z = exp(-NI ) 



dy 



exp 



;4V + ' 



(3) 



Vy 3 + o, 

A r 3! VA^ 



(8) 



The rest of the work is to perform a perturbation with respect to l/N which is also 
called the loop expansion parameter. The loop expansion is well known as the semi- 

(2) 

classical or WKB approximation. In almost all cases, Jj ' / so that l/N expansion can 
be performed with the aid of the Gaussian integration, however there is an interesting 
situation: 1^ = , when A 2 < and |A| = uj/2 , called caustic Utilizing a 

standard prescription with the Airy function in this region and l/N expansion under the 
Gaussian integration in the other region, we can conclude that AFM does work satisfactory 
even within the next leading approximation but more excellently if higher order effects 
would be taken into account. 

In this paper, we pursue the study for checking validity of AFM: models are 1- 
dimensional four-fermi ones with N species, that is, quantum mechanical four-fermi mod- 
els which can be solved analytically. We study two types of model: one has the anti-normal 
ordered form; since whose Hamiltonian in path integral formulas is expressed as the func- 
tion of the Grassmann numbers , not of (the normal ordered form) or of 
+ £?-l)/ 2 ( tri e Weyl ordered form) [TH]. We establish the result that AFM works 
well by taking higher orders in the loop expansion as was the cases in the bosonic 8 and 
the O-dimensional fermi cases jH], which is the content of the section 2. The second model 
is a simpler model of the number operator. Classically there is no difference between these 
two models but here arises an interesting situation: all the higher order corrections seems 
to vanish (although we have checked it up to the two- loop order), which reminds us of 
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the WKB exact models (see the papers ^1] and the references therein). The result that 
the lowest order approximation (because all the higher orders do not contribute) well fits 
to the exact value is also found, which is the content of the section 3. The final section 
is devoted to the discussion where the failure of the path integral representation in the 
continuous time and a new trial toward chiral symmetry breaking in QED are presented. 
In the appendix, the exact energy eigenvalue of the ground state is discussed. 



2 The Model-(l) 

The starting Hamiltonian is 



A 2 

H{a ] , a) = -uj& • a f + — ^ o^ataf , (9) 



a • a 1 " = > ( 10 ) 



where 

A 

((,("■' 

1=1 

a\ and hi (i = 1, 2, • • ■ , N) are the creation and the annihilation operator of i-th fermion, 
satisfying 

{fij, at} = <% , {aj,%} = {a|, a]} = ; (i,j = l,2,---,N). (11) 
Introduce the number operator 



n = a 



N 

t.aU^afol , (12) 

whose eigenstates are 

h\n,r) =n\n,r) ; (n = 0, 1,2, . . . ,N ; r = 1, 2, . . . ,( ^ \\ (13) 

=«l i al 2 ---alj°) ! (*i >*2 >•••>**) ; (14) 

where r specifies the number of combinations of N elements taken n at a time without 
repetition. 

The Hamiltonian @ reads in terms of h as 

*(**,*> = _( w -* + *)* + („-a» + *)» + *1a., (15) 

whose eigenvalue, H(a*, &)\n, r) = E n \n,r) , is 

E n = -(oj- — + —) N+(lu-\ 2 + —) n + — n 2 . (16) 
V 2 2JV7 V 2N J 2N v; 



: We employ the anti-normal ordering in the Hamiltonian @ because in the path integral representation we 
can have matched subscripts between £j and £*: see p|] an d Eq. (J2UJ. 
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The ground state energy Eq is given as the lowest energy state out of N + 1 eigenenergies; 



E G = mmE n , (n = 0, 1, 2, . . . , N) ; 

whose explicit calculations are relegated to the appendix. 

Meanwhile Eq can be picked up from the partition function Z(T), 



(17) 



N 



Z(T) = Tr (e~ TH ) = ^ £ <n,r|e" 



n=0 r=l 



such that 



ffe = - lim -lnZ(T) 



lim — In ( 1^1 

\n=0 v 7 



,—TEn 



(18) 



(19) 



whose final expression apparently coincides with the definition Q17|) . Z(T) has a path 
integral representation, 



Z(T) 



lm vn/7 

£{(S • ^) + At ■ co + ^(5 • ^) 2 }} 



x cxp 



, (20) 



AP 



where Z\i = T/N-t and AP stands for the anti-periodic boundary condition, £o = ~£,N t 
|lUj . Estimating Z(T) under AFM and comparing the results with the exact value, we 
can check validity of AFM. 

Introducing the auxiliary field cTj in terms of the Gaussian identity 



1 = lim TT / \ — da 



'At 



2vr 



Tj CXp 



At 



iX 



(21) 



to erase the (£* -£) 2 term in Eq. Q2U|) and performing the Grassmann integration d N £jd N £j 
we obtain 



Z(T) 



N t 

lim I I 

AT t ->oo - LJ - 



2vr 



dcr,- exp[— iVJ[cr]] , 



where 



J[cr] = ^ of - In det {<% - <5;_i j + At (oj + iXai) <%} . 



(22) 



(23) 



As was explained in the introduction 1/N is the loop expansion parameter. When N goes 
larger (although we will assume that N is not large in the following), it is expected that 
the integral is dominated by the saddle point aoj obeying the equation of motion, 

5I[a] 



5a * 



At(aoj — iXSj 



33) ' 



(24) 



where the Fermion propagator Sij obeys 

N t 



^ {&ik - k-\,k + At (lj + i\(J i) 5 ik } S kj = Si 



(25) 



k=l 



Expand I [a] around <7oj such that 

1 r(2), 



I[a] = I[<t ] + - iy (ai - cr 0i )(fTi - cr 0j ) 



oo j 2V< 
fc=3 ' 31,32,- ,ifc=l 

where we have used the abbreviations, 

<5 fc lM 



Ah) 

31— J'k 



(27) 



<T=(70 



Shifting and scaling the integration variables, we obtain 

Nt . j-TT r n AT f 

Z(T) _ e-^W^n/V^exp -iE^ 



oo 



oo Nt 

E l \ ^ r (k) 

J,|7\rfc/2-l 2^ 1 hh-3k a ii CF - 



k\N k / 2 ~ 

k=3 3i,n,- ,3k=l 



32 ' " a jk 



where 



iff = At A- 1 = At(5ij - X^AQSySfi) 



n 13 

with Aij being a propagator of the auxiliary field, and 

r(3) _ ;\3fA+\3 



iX (At) (SijSjfcSki + SikS^jSji) , 

A (At) 4: (SijSjkSkiSii + SijSjiSikSki + SikSkjSjiSu 
~\~SikSkiSijSji + SuSijSjkSki + SuSikSkjSji) 

Armed with these we have 



ijk 
'ijfcZ 



^2-ioop = exp 



-iVJH-ilndetA^ 1 



1 + — (two-loop graphs) 



(28) 

(29) 

(30) 
(31) 

(32) 



where the two- loop graphs are shown in Fig. ^ the tree part and the tree with the one-loop 
contribution are given as 



y 

^trcc 
■^1-loop 



-NI[a ] . 



exp 



.jV/H-i In det^r. 1 



(33) 
(34) 



respectively. 
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Figure 1: Two-loop graphs: the vertices, I® flHOj) and 1^ ()<U|) are contracting to the points in 
the left hand side to recognize the two-loop explicitly. A, the dotted line, denotes the propagator 
of the auxiliary field. In the right hand side, the line designates the fermion propagator ()25p 
and the numbers upper in (a) ~ (e) denote multiplicity. 



Now let us start a detailed estimation: take a time-independent solution, 

a j i-> a Q ; (35) 

with the overlined symbol, since we are interested in the ground state (=the vacuum). The 
Fermion propagator (|25|) can be calculated with the aid of the anti-periodic eigenfunctions 
/ p (0) = —fp(Nt): 

f p (j) = lj^P+W ; (l<j<N t , l<p<N t ), (36) 



obeying 



as 



N t N t 

E fp(i)f* P ( k ) = 5 3* , E f* P u)m = 6 Pq , (37) 
p=i j=i 



N * i N * ( r \k-j 

s jk = g^(p)/pw=-^g gp _ ( ( r + ^) 

= (i + nit)* + 1 l M1 + f2At ^ Nt+k ' 3 ' 1 ~ W(i + fl^)*-'"- 1 } , (38) 
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where 



and 



Q = U + iXao , 

1 : i>j 
: % < j 

, _ i7r(2p+l)/JVi 



(39) 
(40) 

(41) 



is the N t th root of (—1). The equation of motion (|24|). now called the gap equation, 
becomes with the help of the propagator (|38|) 



<To = iXSu = iX 



In view of Eq. (|3*9*|) , this reads 



[l + DAt)"*- 1 Nt ->oc e OT 
-> iA- 



(1 + fi^t)^* + 1 



(42) 



f2-co = -y 



whose right hand side becomes a step function as T« oo (See Fig. |2J). 



(43) 



















T — > Large 







Figure 2: The shapes of the right hand side of the gap equation ()43j) when A 2 > 0. Those 
to the step function when T i— ► oo. 



The overlined propagator Z\y of the auxiliary field (f2*9*|) can be expressed as 



Tb 



a a(a + Tb) 



fipNt 



(44) 



where 



a = 1 - Z\tA 



2 (l + Jldt) 



N t -2 



(1 + fl^t)^* + 1 



b = A 



2 (1 + ftAfc) 



AT t -2 



{(1 + QAt) N * + l} 2 ' 



(45) 
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and use has been made of the periodic eigenfunctions F p (0) = F p (Nt): 

1 



F P (j) 



-P-jripj/Nt . 



(l<j<N t , l< P <N t ) 



with 



Aft N t 

Y,F p (j)F;(k) = S jk , ^F;(j)F q (j) = 5 pq . 
p=l j=l 

Therefore the tree part is found as 

Z(T) tree = exp(-NI[a ]) ; 

where 



-Oq — lndet{<5y — Si-ij + Atil} 



T 

2A 2 



(f2 - uj) 2 - ln(l + e 



We need to know det A^- in order to obtain the one- loop part (|34j) : 



det A: 



'.J 



1 - AtX 



2 (l + QAt) 



N t -2 



(1 + f2At) N * + 1 



N t -1 



(l + SlAt) N > + l {(l + ft4i) w + l} : 



Nt^oo 



exp 



-A 



e OT + 1 



1 + A 



1)' 



Note the term proportional to At in the second line vanishes when Nt i—* oo. 

The two-loop parts shown in Fig. ^ can be calculated by putting the overlined 
titles, Aij and S tj into ig, (JH0|, and J^, (EU): 



Fig. da) 



A 4 

^(^*) 2 SijSjkSkiSuAitAji 

i,j,k,l 



N t ^oo 



A 4 T 2 e 2i7T 



+ 



A 6 T 3 e 3OT 



+ 



4(e fiT + l) 4 2(e OT + l) 6 

A 8 T 4 (e 2OT _ 4e OT + 1)e 3fiT 



1 + 



1 + 



A 2 re^ T \ 
A 2 Te^ T ' 



Fig.m(b) 



A7 



N t — >oo 



A 4 - - 

--^-(^i) 2 ^ SijSjkSuSuAijA 

i,j,k,l 

A 4 T 2 e 3 ^ T A 6 T 3 (e^ T - l)e 3 ^ T 



A 8 T 4 (e 2fiT _ 4e OT + 1)e 3fiT 



1 + 



A 2 re^ T 

( e l7T + !)2 
-2 



-1 



+ 



12(e^ T + 1) 



1 + 



A 2 Te^ T 
' e DT + 1)2 



A 6 

Fig. Hie) = — — (At) 3 SijSjkSkiSi'jiSjik'Sk'i'AijA kk iAiiji 



Nt^oo 



2 

i,j,k,i' ,j' ,k' 

A 6 T 3 e 4i7T ^ A 8 T 4 (2e^ r - l)e 4 ^ r / A 2 Te^ T 



2( e «r + X )6 2 (e OT + l) 8 ^ (e OT + l) 2 / 

A io r 5 (e flr_ 1)(5e flr_ 1)e 4flr / A 2 Te i7T y 2 
8(e^ T + l) 10 \ + (e° T + l) 2 J 

A 12 T 6 (e i7r_ 1)2e 5J?T / A 2 Te ^T \ ~ 3 

+ 8(e^ T + l) 12 l 1+ (e^ T + l) 2 J ' (53) 



A 6 

Fig.HJd) — — (At) 3 SijSjkSkiSiij>Sjik'Sk>i>AuiAj k iA k ji 



i,j,k,i',j',k' 



A 6^3 e 3i7T 

A 8 r 4 (e 2^T _ 4e J7T + 1)e 2^T / A 2 Te^ T V' 

+ l2(e^ T + l) 8 + (e^ T + l) 2 J 

A 10 T 5 (e OT_ 1)2e 4flr / A 2 Te ^T \- 2 

8 ( e l7T + 1 )10 ^ + ( e OT + !)2 J 

A 12 T 6 (e^ r -l) 2 e 5 ^ / X 2 Te^ T Y* 

24(e^ T + l) 12 I + (e^ T + l) 2 J ( } 



A 6 

Fig. |TJe) = — — (Z\i) 3 SijSjkSkiSi/jiSj/k'Sk'i'AiiiAjj'A kk i 



N t — >oo 



6 

i,j,k,i' ,j' ,k' 
A 6 T 3 (e 2OT + 1)e 2OT 



6 



12(e^ T + 1) 
+ ^ 1 + 



2 



6(e OT + l) 8 \ (e OT + 1) 

A 10 T 5 (e^ T - l) 2 e 4 ^ r / ^ A 2 re^ T \ ~" 

8(e° T + l) 10 \ + (e^ T + l) 2 J 

A 12 T 6 (e^ T - l) 2 e 5 ^ r / , A 2 Te^ r Y* 

+ 24(e^ T + 1) 12 ^ + (e^T + 1)2 J " ( ^ 

(Note that a difference between Fig.U(d) and (e); the propagator Aj k iA k y in (d) becomes 
Ajj>A kk > in (e).) 
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According to the definition ()19|) we obtain 



j-i2-loop 



N 



as well as 



pitree 



E, 



1-loop 
G 



lim — In ( 1 + e 



lim 



A 2 e fiT 



+ lim 



1 



T^oo 2 e ftT + 1 T^oo IT 



-— lim 

iV T^oo 



A 4 T(2e m - l)e : 
4(e^ T + l) 4 



J2QT 



A b T 2 (7e 



He«r + 1)e 2fir 



+- 



12(e OT + l) 6 

5A 8 T 3 (e fiT _ 1)2e 3OT 



N 
N 



24(e OT + l) s 



1 + 




lim — In ( 1 + e 



SIT 



sit 



lim 



A 2 e nT 



1 -^oo 2 e fiT + 1 T^oo 2T 



lim — In ( 1 + e 

T->oo T 



1 , / X 2 Te^ T 
+ nm — — m 1 + 



+ 1) 



(56) 



(57) 



(58) 



Estimation is made for two cases that (i) A 2 > and (ii) A 2 < 0. 

(i) A 2 > : the solution of the gap equation Q43JI can be categorized into three cases, (ia), 
(ib), and (ic), according to the value of lo: see Fig. |31 



Figure 3: The gap equation: A 2 > : the step function stands for the right hand side of the 
gap equation ()43j) and the line, Q — oj , for the left hand side. There is one crossing point in 
each case: (ia) ; oj > and A 2 > oj . (ib) ; oj > and < A 2 < oj . (ic) ; oj < . 



(ia) to > and A 2 > to : ft = then from Eq. lj4*3|> . 



QT 



LO 



A 2 
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^2-loop = _ N _^ ,S i A L 



(59) 



so that from Eq. (|56|) 

2 A 2 2 JV 

whose O(N), O(l), and 0(1/N) terms designate the tree, one-loop, and two-loop 
contributions, respectively. 

(ib) u) > and < A 2 < u> : J? = u; - A 2 , then 

e 3 " ► OO , 

so that there is no correction from the two-loop in Eq. to give 



E, 



l-loop 
G 



E, 



2-loop 
G 



(60) 



whose 0{N) and 0(1) terms are the tree and the one-loop contributions. There are 
two domains in A 2 divided by u>, which exposes a striking difference to the exact 
case (|1 15|> where N + 1 domains exist. Due to these, the curve of exact energy has 
many discontinuities arising from the boundaries of domains. We plot the results 
of (ia) and (ib) with N = 2, u = 1 in the left of Fig. |1J Although the deviations 
in the tree and one-loop results from the exact energy is substantial, the two-loop 
contribution cures the situation; whose effect is much clearly seen with N = 4 in the 
right of Fig. EJ 



0.5 

-0 . 5 
-1 

-1.5 
-2 

-2.5 




Figure 4: The case of A 2 > ; u = 1, corresponding to (ia) and (ib). The solid, dotted, 
dashed-dotted, and dashed-double dotted lines designate the exact, tree, one-loop, and two- 
loop results, respectively. We put iV = 2 (left), 4 (right). The approximation becomes better 
when taking higher loops into account. 

(ic) oj < : Q = — \u\ , then 

so that the one- and two-loop results vanish, yielding to 

E% ee = ^- loop = E 2 G loop = , (61) 
whose result coincides with the exact energy (|117|) in the appendix. 
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Figure 5: The gap equation, A 2 < : (iia) ; uj > . (iib) ; uj < and < |A 2 | < \u\ . ( 
uj < and |A 2 | > \uj\ , where three solutions, (iicl) ~ (iic3) , exist. 

(ii) A 2 < : five different solutions of Eq. (J1H|) are found according the value of uj. See 
Fig. El 



(iia) uj > : Q = uj + \ \ 2 \ , then 



e — > oo 



so that there is also no correction from the two-loop, yielding to 

U-loop _ T-.2-100P / |A | \ | A 



s ™ = E ^o P = _ N | u + r_> + r_, > (62) 



which is nothing but the result obtained in Eq. 1)118(1 in the appendix. The tree 
result has a slight deviation but the one-loop correction fits in the curve. (See the 
left graph in Fig. El) 

(iib) uj < and < |A 2 | < \uj\ : Q = -\uj\ , then 

so that in Eq. (jSHj) the one- and two-loop corrections vanish to give 

£g ee = £g 10 ° P = ^G l0 ° P = • ( 63 ) 

(iic) uj < and |A 2 | > \uj\ : from Fig. El there are three different solutions, 
(iicl) Q = — \uj\, then 

so that there is again no correction from the higher orders to obtain 

£g ee = £g 10 ° P = ^G l0 ° P = • ( 64 ) 

(iic2) Q = , then from Eq. 



|A 2 I - U 
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so that 

^2-loOp »r |^| 2 . M . 1 |^| ,,, r . 

^ =iV ^ + T + iv— ' (65) 



from Eq. (|56j). which is apparently positive to be greater than (iicl) and should 
be discarded. 
(iic3) Q = -\uj\ + |A 2 |(> 0) , then 



QT T-*oo 
e — ► oo ; 

so that in Eq. (|56j) the two- loop correction vanishes, yielding to 

ix2i 2iV 1 

positive : < A < M 

N — 1 

,x2, 2iV , 

negative : A > \uj\ 

11 N _ 1 i 1 



(66) 




Figure 6: Left: A 2 < ; u; = 1 case, corresponding to (iia). Here also N = 2. The solid, 
dotted, dashed-dotted, dashed-double dotted lines designate the exact, tree, one-loop, and two- 
loop results, respectively. The tree approximation is slightly deviating from the exact value 
but including the higher orders reproduces it: the one-loop approximation is accurate enough. 
Right: A 2 < ; = -1 , and N = 2 case, corresponding to (iib) and (iic) : the two- as well as 
the one-loop approximation matches with the exact value. 

2N 

Therefore we should adopt the solution (iicl) up to |A 2 | = — — -|u;| then switch it to 

(iic3), whose result is plotted in the right graph of Fig. with putting N = 2, u = — 1. 
After including the one-loop effect, the ground state energy as well as the number of 
domains coincides with the exact energy ()119j) in the appendix. 

So far the tree or one-loop result matches with the exact curve except the case A 2 > 
; uj = 1 (ia) and (ib), where the number of domains (two) are different from the exact one 
(N + 1), Eq. (jll5|) . The two-loop effect, however, cures the situation and the deviation 
becomes smaller in N = 4 (Fig. |3J). The reason is seen from Fig. [7| where we can recognize 
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that the tree approximation approaches closer to the exact value when N moves from 2 
to 4 and 10, and moreover that when TV" goes larger, discontinuities in the exact curve 
fade away. (As the number of domains in A 2 Q115JI increases sharpness at the boundary 
wears off gradually.) Needless to say that AFM is analytic so that unless a boundary of 
domains in the loop expansion happens to coincide as of A 2 < ; u < 0, (iib) and (iic) 
(Fig. EJright)), the deviation is inevitable as of A 2 > ; to > 0, (ia) and (ib). Therefore 
we conclude that AFM works still well even in A 2 > ; u> > : the deviation emerges not 
from the failure of AFM but from the model with discontinuities in the energy curve. 



exact - 
tree - 



N = 2 



N = 10 



exact - 
tree - 




Figure 7: The exact (the solid line) and tree (the dotted line) results for iV = 2, 4 and 10 , by 
putting A 2 > ; uj = 1. Discontinuities in the exact energy curve wear off gradually. Although 
the vertical scale (the energy value) is different, the tree results come closer and closer to the 
exact value. 



3 The Model-(2) 

In this section, we adopt a slightly different model: 

A 2 A 2 
H{a\a) = u/(fit . a) + — {a) • o)(o t • d) = un + — h 2 . (67) 

Classically, there is no difference between the model-(l) and -(2), but the reason for 
considering this model is that all the higher order corrections seem to vanish (although 
we have checked up to the twodoop) to give us another example of WKB exact model |14j . 
The energy eigenvalue is obtained as before, 

H(a\a)\n,r) = E n \n,r) ; 
A 2 

E n =iun+ —n 2 . (68) 
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The minimum (= the ground state) energy Q17|) . calculated in the appendix, should be 
compared with that obtained from AFM: the partition function is 



r r r ( 



(69) 



AP 



which, by use of the auxiliary field, turns into a new form , 

JV, 

N t -*oo 



Z(T) 



lim 1 I 



2tt 



d<Tj exp[— iVJ[cr]] , 



(70) 



where 



I[a) 



T[u + 



+ ~Y ~~ lndet ~~ + <^ ^ w + A 2 + *M - T^y^ j 



(71) 



The saddle point doj fulfills the equation of motion 



Z\t(<7oj — iXSjj) 



(72) 



f=CT0 



where obeys 



J — "ij ■ 



(73) 



Note that contrary to the previous situation, 0(1/N) term has already appeared in the 
expression. Therefore the loop expansion should be distinguished from the 1/N expansion 
in this model: we should perform the loop expansion then arrange the results in the order 
of 1/N. All the procedures to the 0(1/N) (two-loop) ground state energy are, however, 
the same except that all quantities are tilded; the time-independent solution of the gap 
equation (|4*3~)) should be expressed such that aoi t— > d$ . 

Therefore the Fermion propagator, the solution to Eq. I|73|) . is given as 



Sjh = s-^ {0 jk (l + QAt) 

(1 + DAt)Nt + l I 



where 



Nt+k-j-l 



A 2 



1 Ki 



^(l + nAt)*-*- 1 } , (74) 



fl = lo + A + i\a - — 



(75) 
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With the aid of the propagator (|74|). the equation of motion 1)72(1 yields to the gap equation, 

do = lADa = iX = — > iX— . (7b) 

(l + nAt) N * + l e aT + l 



Namely, 



A 2 



e 



QT 



Q-[uj + X I - — )=-X z ^ ; (77) 

V 2N) e a T + 1 > V > 

m view of Eq. lf75|). In terms of l/N expansion, Q can be expressed as 

^ = ^o + f + % + • (78) 

Inserting this into the gap equation (|77|) and arranging the power of l/N, we obtain 

p A)T 

^o-(- + A 2 ) = -A 2 ^ TT , (79) 

and 

\2 / p D T X" 1 

ft = ~[l + TX \ ^ , - 1 , (80) 



2 I (e fi o T + r 



2 y ( e ^oT + !)3 ^ ( e ^ r + X )2 

The tree part of the partition function (|33|) is 

Z(T) tree = exp(-JVI[fr ]) ; (82) 

with 

f X 2 \ T - 
I[a ] = Tlu + — )+-d%-lndet{5 ij -5 i - 1J +Atfi} 

\2\ T / _ \2\ 2 



= T ^ + Tj-2A5(«- W - A + m) -Ml + e'"). (83) 

The one- loop partition function (|34j) with Eq. (|5Ujl and the two- loop one shown in Figs.^ 
(a)-(e) or Eqs. (|51|) ~ (|55|1. are all the same with the substitution I? i— ► I?. 

The ground state energy up to the two-loop is obtained from the second to the fourth 
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line in Eq. (|56|) as well as from Eq. (|83|) such that 



E, 



2-loop 



N 



lim — ln(l + e 



A 2 e QT 1 / 

lim = h lim — In 1 + 

T^oo 2 e l?T _(_ i T-*oo 2T \ 



— lim 

N T^oo 



A 4 T(2e^ T - l)e : 
4 ( e 4r + !)4 



,2I2T 



1 + 



A 6 T 2 (7( 



,2J?T 



lle OT + l)e 



2QT 



12(e^ T + l) 6 



1 + 



A 2 Te^ T 
(e^ T + l) 2 

( e OT + X)2 J 

X 2 Te^ T 



-i 



+ 



5A 8 T 3 (< 



SIT 



l) 2 e : 



2 Q 3OT 



24( e «T + X )s 



1 + 



\ 2 Te^ T 
(e^ T + l) 2 , 



(e^ T + l) 2 , 

-3 



Expand J? in terms of 1/N according to Eq. ()78|) to give 



1/jV 



N 



- + |-^»-(- + A 2 )) 2 - T l-iln(l 



+ lim 



1 



T^oo 2T 



— lim 



In 1 + 



A 2 Te^° T 
>«oT + 1)2 



A 4 T(e 2 ^ T - 4e D ° T + l)e^ T / A 2 Te^ T 

\ + (e n ° T + l) 2 

A 2 Te^oT V 2 

1 + 



e n T + X )4 

A 6 r 2 (4e 2i7 T _ lle i7 T + 4 ) e 2A,r 
12(eA,T + !)6 

5A 8 T 3 (e^ T - l) 2 e 3 ^ r 
24(e^oT + i)8 



1 + 



A 2 re^° T 
( e «oT + !)2 



( e a,T + !) 2 

-3 . 



where the first line is the leading approximation, 



Eg = N 



0J + 



A 2 



1 



2\\2 



{Q Q - {u> + A )) 



2 2A 2 

and the term including 0(1) approximation is 



lim — In ( 1 + e 



,o t 



E, 



G 



N 



. + f-^(«o-(. + A^- T l^iln( 1+ e*- 



+ lim 



1 



T^oo 2T 



In 1 + 



A 2 Te^ T 



Again estimation is made for two cases that (i) A 2 > and (ii) A 2 < 0. 
(i) A 2 > : the solution of the gap equation ()77j) can be categorized into three cases, 
(ib), and (ic), according to the value of ui: see the left of Fig. [HI 
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Figure 8: Left: the gap equation ; A 2 > 0: there is one crossing point in each case: (ia) ; u> > 0. 
(ib) ; u < and A 2 > \u\. (ic) ; uj < and < A 2 < \uj\. Right: the ground state energy 

for N = 2: A 2 > ; u — — 1, corresponding to (ib) and (ic). All the corrections vanish and the 
lowest result (the dotted line) almost reproduces the exact value (the solid line). 



(ia) u) > : = uj then from Eq. (f?9"jl 

so that all higher orders of 1/N in Eq. (|85[) do vanish, leaving the leading O(N) 
term, 

^G - ^G - h G - U ' 
which matches with the exact energy (|122|) in the appendix. 

(ib) u < and A 2 > |w| : J? = then from Eq. $7^ 

| a; | 

so that there remains only the leading term in Eq. (|85[) to give 

(ic) oj < and < A 2 < case: J? = — M + A 2 (< 0) then 

e n T T->co Q ^ 

so that the leading term only gives us 

Eg = Eh = E 1 c ( N = n(-\u\ + ¥-) . (90) 



2 , 

We plot the result in the right of Fig.|S]from which we see that the leading approximation 
almost matches the exact curve. The deviation comes from the boundary of A 2 domains, 
where discontinuity in the exact energy curve emerges. Note that the number of domains 
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Figure 9: The gap equation, A 2 < 0: (iia) ; u > and < |A 2 | < uj. (iib) ; u > and |A 2 | > 
where three solutions, (iibl) ~ (iib3) , exist, (iic) ; u < 0. 

is three in the energy eigenvalue (generally N +1: see Eq. (|127|> in the appendix) but two 
in AFM. 

(ii) A 2 < 0. Five different solutions of Eq. (|T7|) are found according the value of uj. See 



Fig. El 

(iia) uj > and < |A 2 



< ui : i?o = w then 




oo 



so that the leading term in Eq. I)85|) only gives us 



TP" — TP 1 — TP 

&cx — &a — tin 



,1/JV 
G 



= . 



(91) 



(iib) uj > and |A 2 | > u : 
(iibl) i?o = uj then 




oo 



so that again the leading term gives us 



Eq = Eq = E G /N = . 



(92) 



(iib2) Qq = then from Eq. (S) 




so that there is no correction other than the leading term yielding to 



,1/N 
G 



,2 



TpN _ TP 1 — TP 
hin — tLn — Hi. 



= N 



2IA 2 



(93) 



which is positive to be greater than (iibl) and should be discarded. 



20 



(iib3) f2 = co - |A 2 |(< 0) then 



e ti T T-+CO Q ^ 



so that there remains only the leading term in Eq. (|85[) to give 

at 1 i/)v / IA 2 |\ f positive: < |A 2 | < 2uj 
Eg = E^ = E l J N = N [u - i^i ) = { ' . (94) 

V 2 / [ negative : |A 2 | > 2u 

Therefore we should adopt the solution (iia) up to |A 2 | = 2uj then switch it to (iib3): 



Eg 



: |A 2 |>2w 



2 / , (95) 



: < |A 2 | < 2uj 



which is nothing but the exact energy (|129|) in the appendix, 
(he) u> < : Qq = — \uj\ — |A 2 | then 

e«» TT =?0, 

so that there is no correction from the higher orders to obtain 

I A 2 



E^ = Eh = E 1 c / N = -N[\u ) \ + ^) , (96) 



2 

matching with the exact energy (|13U|) in the appendix. 

We have checked the model up to the 0(l/iV 2 ) to find that all higher order corrections 
vanish as well as that the leading term reproduces the exact value except the case A 2 >0 
; u < 0. Therefore we can conclude that the model belongs to the WKB exact class |14j . 
although the expansion is performed with respect to 1/N which differs from the loop 
expansion(= WKB) in this case. A slight discrepancy emerges from the boundaries where 
discontinuity is eminent (The right of Fig. 0). As was stated before AFM is analytic 
so that if the number of the domains in A 2 differs each other a slight deviation could 
inevitably occur. 



4 Discussion 

In this paper, we discuss validity of AFM in terms of quantum mechanical four-fermi 
models. The model with the anti- normal ordered form shows us that when A 2 > ; uj > 
the two- loop result almost fits the exact value (Fig. HJ, when A 2 > ; uj < the 
tree result becomes exact, and when A 2 < ; u > (Fig. EJleft)) as well as when 
A 2 < ; u < (Fig. Efright)) the one-loop correction reproduces the exact energy. Even 
m A 2 > ; uj > 0, we should expect the exact fit of the one- loop or tree result, but the 
non-analytic structure of the exact energy curve causes the deviation. 

In the second model, all the higher order corrections vanish other than the lowest one, 
showing us another example of the WKB exact class: when A 2 > ; uj > the leading 
term reproduces the exact value, when A 2 > ; uj < it fits well except the region 
where non-analytic structure becomes eminent (Fig. fright)), and when A 2 < ; uj > 
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(Eq. (|94|)) as well as when A 2 < ; to < (Eq. (|96j0 it fits exactly. The only deviation in 
(Fig. IBtright)) comes from the non-analytic structure of the exact energy curve. 



So far discussions are made under the discrete time path integral representation which 
keep Nt finite until the end of calculations, but in perturbation or WKB approximation 
people rely on the continuous time path integral which takes i — >• oo first. Of course, 
it is simpler and easier, but in some cases |15j we are forced to use the discrete time 
representation. Here we give another example: the model-(l) with Nt i— > oo. The partition 
function (1201) becomes 



Z(T) 



V N £V N C exp 



T 



(97) 



AP 



where AP implies £(0) 



-£(T) . After introducing the auxiliary field a(t) it gives 
Z(T) = J T>aexp[-NI[a\] , (98) 



where 







dt—^-Trlni — +uj + i\a(t) ) , 



(99) 



with Tr being the functional trace. Find the constant solution a(t) = ctq, yielding the gap 
equation 



6I[a) 



8a(t) 



o"o — i\S(t, t) = , 



and the Fermion propagator, 

1 



5(t,0 



2 cosh Of L 



with 



Q = to + iXctq . 
By noting 0(0) = 1/2, Eq. (tTUUl) reads 



A 2 QT 

— tanh 

2 2 



Q — uj 

The propagator of the auxiliary field Q44|) turns out to be 

A( P ) = l 



8 p0 , p=—- (m = 0,±1,±2,---) ; 



1 + Tb 



(100) 



(101) 



(102) 



(103) 



(104) 
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with 



, A 2 / lS QT 

b = — 1 - tanh 2 

4 2 



(105) 



The ground state energy up to the two-loop ( Eq. (|56|) under the discrete time ) is therefore 
given as 



E, 



i2-loop 
G 



N 1 QT 

. — {fl-uif-N^-h,^ — 



+ — lim —In 

2 T^oo T 



— — lim 



TA 2 / , 9 QT\ 

— l-tanh 2 — l-3tanh 2 — 



A 2 T / 2 i7T 

x H 1 - tanh 2 

1 4 V 2 



+- 



A 6 r 2 

768 



1 - tanh" 



QT 



3 - 19tanb 

2 



9 OT\ 

— 

2 y 



A 2 T / , 2 

x H 1 - tanh 2 

4 V 2 



H tanh 2 1 - tanh 2 

1536 2 V 2 / 



x 1 + 



A 2 T 



1 - tanh 2 



QT 



4 V 2 " • (1 ° 6) 

We estimate the case A 2 > and u > 0. The solution of the gap equation (|lU3j) can be 





A 2 

2 













A- 




2 


/SIS) 



exact 
tree 



JV = 2 



Figure 10: Left: the gap equation in continuum. Right: the ground state energy for TV = 2 
(A 2 > ; uj = 1). The exact curve and the approximation result (only the lowest is left intact) 
are given by the solid and the dotted lines, respectively. Disagreement is prominent. 

categorized into two cases, (i) and (ii), according to the value of uj: see Fig. ITUl 

A 2 

(i) uj > and < A 2 < 2lu : Q = oj - —, then 

QT T^oo 

tanh — ^— — ► —1 , 
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so that 



^tree = ^1-loop = ^2-loop = N - ^ (107) 

(ii) to > and A 2 > 2to : I? = 0, then 

. QT t->oo 2uj 
tanh— -> --, 

so that 

£g ce = ^- 1o °p = E 2 G loop = (108) 

In continuum all the higher order corrections vanish. The result for N = 2 and uo = 1 is 
plotted in Fig. 1101 disagreement with the exact value is apparent. In the discrete time (in 
Fig. the tree result is improved by the higher loop effects but there is no higher loops. 

There is a long history of studies in chiral symmetry breaking in QED |16j : in a strong 
coupling region there seems to exist a chiral breaking phase. However trials, such as 
the Schwinger-Dyson or the effective potential approach, have been suffered from gauge 
dependence. Here is an alternative: a starting Lagrangian of massless QED reads, 

C = -^7 M (dp - ieA^) iP - ; = d^A v - d v A^ ; (109) 

so that the partition function, 



Z = J P^Whp U d 4 xCj ; (110) 
can be, with the aid of the auxiliary fields, reduced to 

Z = J VB^VC^VtpVipexp (J d 4 x£'^j ; (111) 
C! = —ihn {dp - ieBf, + e^C^) ip - ^— ^ | (ipip) - (ip^ip) ) 

~~ rM -n + ™ 2 J B »--2 G i>> ( D = d ^ (112) 

where use has been made of the Fierz identity. Note that there is no gauge fixing, the 
remnant of which can be seen from the invariance, \— > B^ + dnA. The four-fermi form 
looks similar to the Nambu-Jona-Lasinio model so that the chiral structure that whether 
the quantity, {ipip), is zero or nonzero would be explored in a parallel manner. This is our 
next task. 



This work is supported by Grant-in- Aid for Science Research from the Japan Ministry of 
Education, Science and Culture; 13135217. 
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A The exact ground state energy 

In this appendix calculations are made for the ground state energy (|17|). 



A.l Model-(l) 

First complete the square in the energy eigenvalue (|lfi|) to give 

A 2 f / 1 Nuj\) 2 N ( A 2 \ . . 

The vertex is thus 

1 Nu . . 

n vert = N - - - , (114) 

where the minimum or maximum occurs if there is no restriction for n. However < n < 
N so we need a detailed inspection. 

(i) A 2 > and to > : E n is concave and n vert < . The pattern of E n is given 



in Fig. El From Fig. HTf a). if 



1 / N 



n V ert <-<^0<X<uj 



2 \N-1 
n = is the minimum to give 

E G = E = ^Tp^A 2 - Nlo . 

Generally, from Fig. 111( c) , if 

2m - 1 2m + 1 / N \ ( N 



< n vert < — — — 44> U) ( — ) < A 2 < UJ 



N — m) \ N — m — 1 

n = m(l, 2, ... N — 1) is the minimum to give 

(N - m)(N - m - 1) >2 , AT 

Cases can be summarized to give 

" f TT~ — ) < A 2 < u ( N ) , (m = 0,l,2,...JV-l) ; (115) 
\JM — m J \N—m — lJ 

then 

(N — m)(N — m — 1) ~ . T , . 

E G = E m = [ - L -A 2 - (AT - m)w . (116) 

(Note that when m = the lower bound u> of A 2 should be discarded.) As 
goes large the number of domains given by Eq. 1)115(1 increases, which smoothes the 
discontinuity of the energy curve around the boundaries as seen from Fig. with 
N = 2,5 and 10. 
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Figure 11: The patterns of the ground state energy for A 2 > 0. (a) shows the case of n vert < -. 
The minimum occurs at n = 0. (b) shows the case of n vert > . The minimum occurs 



at n = N. (c) shows the case of - < n VCTt < — — 
1, 2, . . . , N — 1) the minimum occurs at n = m. 



2 

1 2m - 1 2m + 1 , 

— . When — - — < n vert < — - — ; (m 



(ii) A 2 > and uj < : 



n V ert =N-- + -4^- > N - - . 



2 A 2 ' " 2 
From Fig. fTTTb). n = N is the lowest. E N = from Eq. (fTT3|) . so that 



-Eg = En = . 



;ii7) 



(iii) A 2 < and u > : E n is convex and the pattern of E n is shown in Fig. 1121 

1 Nlo N 

since N > 2. Therefore from Fig. ll2f bh n = is the minimum and Eq = — |A 2 | — Nuj 

so that 



N-l, 



|A 2 | — Nlo . 



(118) 
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(a) 



Figure 12: The patterns of the ground state energy for A 2 < 0. (a) shows the case of n ver t < 

N 

The minimum occurs at n = N. (b): the case of — < n vert - The minimum occurs at n — 0. 



N 



(iv) A 2 < and lo < : 
From Fig.lEta), if 



"vert = N 



1 N\u\ 



2 |A 2 | 



N ,,,, , 

n VPr t < — < A z < a; 



2iV 



2 " " ' '" 1 - , ~ i N- 1 ' 
n = A" is the minimum to give Ejy = 0. Further from Fig. I12lf b) . if 



N 



2N 



— < n vert & \u, N _ _ 



< |A 2 | ; 



n = is the minimum to give £7o 





N - 1 



A 2 1 + TV | a; | . Therefore 



#G = < 



< I A 2 1 - 



2N 
N — 1 



•— 2— |A 2 |+iVM : |A 2 | > ^ryl^l 



A.2 Model-(2) 

Again complete the square in the energy eigenvalue (|68|) to yield 



A 2 ( Nu\ 2 Nlo 2 



E " 2 A I" ' A 2 ; 2A- 



Then the vertex is 



^-vert 



Nlo 



(119) 



(120) 



(121) 
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(i) A 2 > and to > : E n is concave. If 



Nuj 



n ver t = -Tn" < ; 



from Fig. Illf a) n = is the minimum to give 

Eg = E = . 

(ii) A 2 > and uj < : 

^"vert ^2 ' 

which is positive. From Fig. Illf a). if 

1 9 

n ve rt < - 44> A > 2N\u\ 

the minimum occurs at n = to give 

E G = E = . 

Generally if 



2m - 1 



2m + 1 2AT 
< n V ert < z ^ - r M < A < 



2N 



2 ~ 2 2m + 1' 1 ~ 2m- 1' 

from Fig. Hir e), n = m(m = 1, 2, . . . , N — 1) is the minimum to give 

A 2 



E G = E m = -\uj\m + 



When n vert > 



2N- 1 



< A < 



2N 



122) 



(123) 



(124) 



(125) 



(126) 



2N -V 

from Fig. 111( b). the minimum occurs at n = N to give 

iVA 2 

E G = E N = -N\u\ + —- . 



In view of Eqs. (|123|) . (|124f) . and (|126|) . we can take the ground state energy as 
Eq. CEEJ 



A 2 

E G = E m = -\uj\m + —m 2 ; m = 0, 1,...,N 



in the region 



2N 



2m + 1' 



\u\ < A < 



2iV 



2m - 1 



I a; I ; m = 0, 1, . . . , N ; 



(127) 



(128) 



under the condition 



2N 



2m + 1 



m=N 



, 2iV 

and w 

2m - 1 1 1 



oo. 



m=0 
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(iii) A 2 < and to > : E n is convex. 

Nu 

Invert — | ^2 1 ' 

The situation is similar to (iv) in Model- (1). If 

N ,,2, 
™vert < y ^ |A \ >2U 

then from Fig. I12f a). n = N is the minimum to give 

/ |\2| 

E G = E N = N Itu - L-i 



Further if 



y < "-vert ^ < |A Z | < 2uj 



from Fig. 112( b). n = is the minimum to give 

E G = E = . 

Therefore 



Eg 



iVfw- — ) : |A 2 |>2w 



2 7 - . (129) 



: < |A 2 | < 2uj 



(iv) A 2 < and u < : 

^-vert 



iV|w| 



|A 2 | ' 

which is always negative so that from Fig. I12f a) n = N is the minimum: 



L\; Ly - A" ( M + ^ ) ■ (l™) 
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